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The optical properties of the 2x1 reconstruction of the dia- 
mond (111) surface are investigated. The electronic structure 
and optical properties of the surface are studied using a mi- 
croscopic tight-binding approach. We calculate the dielectric 
response describing the surface region and investigate the ori- 
gin of the electronic transitions involving surface and bulk 
states. A large anisotropy in the surface dielectric response 
appears as a consequence of the asymmetric reconstruction 
on the surface plane, which gives rise to the zigzag Pandey 
chains. The results are presented in terms of the reflectance 
anisotropy and electron energy loss spectra. While our results 
are in good agreement with available experimental data, ad- 
ditional experiments are proposed in order to unambiguously 
determine the surface electronic structure of this interesting 
surface. 

PACS numbers: 78.66.Db, 73.20.At, 78.66.-w, 73.20.-r 



I. INTRODUCTION 

Apart from being of fundamental interest, the charac- 
terization of the low-index diamond surface is very im- 
portant from a technological point of view. The fast de- 
velopment of chemical vapor deposition techniques has 
increased the demand for a better understanding of the. 
ground state and excited properties of these surfaces.EJ 
Indeed, much experimental and theoretical attention has 
been paid to the characterization of the geometrical 
structure, vibrational modes and electronic properties of 
these surfaces, with interesting and sometimes contro- 
versial results. In this work, we are interested in char- 
acterizing the optical response of the (111) diamond sur- 
face, and investigating how these properties are related 
to the structural reconstruction and its accompanying 
electronic structure. 

Our interest on this particular 

surface includes concerns on discrepancies between the 
present experimentaMT] and theoretical results.ETEil Ex- 
perimentally, a great deal of the surface electronic struc- 
ture is well known through j-aiigle resolved photoemi&. 
sion spectroscopy (ARPES)pi3 soft-x-rays absorption,l3 
inverse photoemippdpnJ3 and electron energy loss spec- 
troscopy (EELS)D'EI The p^ton-induced measurements 
show a variety of occupiedau and unoccupiecE^Ll surface 
states lying in the fundamental gap. However, the com- 
plete description of these surface states has been difficult. 



since only their dispersion along a few of the main direc- 
tions of the surface unit cell have been measured. On 
the other hand, EELS measurement^ show a prominent 
broad feature at about 2.1 eV, which is attributed to 
transitions from occupied to unoccupied surface states. 
Since EELS experiments measure a transition energy 
which is generally smaller than the difference between 
occupied and unoccupied states, and a relatively large 
uncertainty (±0.6 eV) accompanied this particular EELS 
experiment, a direct comparison with other results has 
not been possible. 

Several theoretical studies have been done to elucidate 
the structural and electronic jpmperties of the C(lllAi 
2x1 surface. Both, ah initim^j and semi-empiricalllil 
theoretical approaches have been employed, yielding 
some differences among them and with experimental re- 
sults. Some of these differences arise from the methods 
ology employed. For example, larlori and co-workerso 
employed a LDA formalism using a plane wave ba- 
sis, while Vanderbilt and Louie, 13 and Alfonso and co- 
workeralj used the LDA formalism based on a set of 
localized orbitals. Lu the former work the energy gaps 
are underestimated, □ as is common in this kind of ap- 
proximation, and a direct comparison to experimental 
results is difficult- JIhe latter theoretical works compare 
well among themja'ty although a systematic shift of about 
1 eV is found when the surface gtfrt^ are compared with 
those measured experimentallycrBflfl On the other hand, 
the semi-ernpirical tight-binding approach of Davidson 
and PickettEll compares well with the ah initio results 
described above, except for an extra shift of the surface 
states by about 0.8 eV. Since the surface states determine 
the location of the Fermi level, there is a large discrep- 
ancy among different theoretical works as to the relative 
position of the Fermi level and the top of the valence 
band that goes from -1.3 to 2 eV. On the other hand, in 
all theoretical and experimental results there is general 
good agreement on the energy gap between empty and 
full surface states. From these considerations, one can 
then conclude that the calculated electronic structures 
alone are not able to uniquely determine the nature of 
the transitions observed in EELS,D'EI resulting in contro- 
versial iiiteijBcetations of the available experimental opti- 
cal data.l3l3'EH3 The evaluation of the surface dielectric re- 
sponse function for this system, and its analysis in terms 
of the associated electronic level structure gives further 
insights into this problem, as we discuss below. 

In the present work, and in close connection with 
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the general description of the optical properties of the 
C(lll)-2xl surface, we investigate in detail the origin 
of the electronic transitions related to the surface re- 
construction. Our calculations employ a semi-empirical 
tight-binding approach that has been used previously 
to study|-the optical ppaperties of Si(lll), (110) and 
(100),mla and C(001)lla surfaces. Our tight-binding 
formalisiu is similar to the one used by Davidson and 
Pickett, Ell except that our extended orbital basis allows 
perhaps a better description of the conduction band due 
to its additional s* orbital, and our level structure is in 
general better agreement with experimental findings and 
other calculations. 

Since we calculate the surface dielectric tensor, 
the results presented here can be compared directly 
with those measured using various optical spectro- 
scopies. In particular, the differential reflectance and 
reflectance anisotropy spectroscopies provide accurate 
information about surface properties of metalsEZI and 
semiconductors Jla This is very important since several 
semiconductor surfaces show a metallic-like behavior due 
to the narrow gap between occupied and unoccupied sur- 
face states that hamper the use of electronic spectro- 
scopies like ARPES, EELS and Scanning Tunnelling Mi- 
croscopy (STM). This indeed seems to be the case for 
the C(lll)-2xl surface, where a narrow gap has been 
found thcoEetiC|ally along one of the main directions on 
the surfacejjilj while no direct experimental evidence is 
found in the literature for this metallic behavior. This 
behavior is in comparison with the Si and Ge (lll)-2xl 
reconstructed surfaces jl3'c2l where the degeneracy of the 
surface states is broken by the buckling of the surface 
atoms, as the theoretical and experimental description of 
the surface states shows. Notice, furthermore, that the 
optical spectroscopies mentioned above have the advan- 
tage over other techniques of allowing in situ real-time 
measurements, which provide the invaluable opportunity 
of monitoring the chemical jSiapor deposition and molec- 
ular beam epitaxial growth£ll as well as the dynamics of 
the chemisorption processE3 The results presented here 
then not only provide answers to fundamental questions, 
but give important information for applications, which 
we expect will motivate future work in this direction. 

In Section II, we present a brief discussion of the struc- 
tural model of the surface and the methods used to calcu- 
late its electronic and optical properties. In Section III, 
we discuss our results and compare with the available 
data in the literature. The results for the optical prop- 
erties are presented in terms of the dielectric response 
of the surface, and the calculated reflectance anisotropy 
and EELS spectra. 

II. MODELS AND METHOD OF CALCULATION 

The diamond (lll)-2xl surface was modeled using a 
slab of 28 C layers with inversion symmetry, yielding a 



free reconstructed surface on each face of the slab. The 
thickness of the slab is large enough to decouple the sur- 
face states at the top and bottom surfaces of the slab. In 
Fig. 0, we show (a) the top view of the surface unit cell 
that contains two C atoms per layer, (b) a side view with 
only the six outermost layers of the slab and (c) the irre- 
ducible surface Brillouin zone (SBZ). Periodic boundary 
conditions were employed parallel to the surface of the 
slab to effectively model a two-dimensional crystal sys- 
tem. The top (and bottom) layer of the slab, shown in 
Fig. |l| withjWack circles, resemble the structure reported 
by Pandey.c2l In this Pandey chain model, the atoms of 
the top layer form a zigzag chain along one of the main 
directions on the surface plane (the x axis in Fig. |l|). 

The coordinates for the six outermost, layers on each 
side were obtained by Alfonso et al.]l3 using a first- 
principles density functional based molecular dynamics 
technique due to Sankey and co-workers .cil (The remain- 
ing central layers have bulk geometry.) The method has 
been employed successfully in studying covalent systems 
such as silicon and carbon. Eil tj The relaxed C(lll)-2xl 
surface obtained with this method showed the zigzag-like 
chains with no buckling on the surface layer, and with 
CC bondlengths of about 1.44 A. The results of Ref. 0] 
are in excellent agreement |J}p±|h previous self-consistent 
first-principles calculations,BtZI where they find unbuck- 
led surface chains with bondlengths equal to 1.47 and 
1.44 A, respectively. The reader is referred to Ref. |2j] 
for a comprehensive description of this technique, and to 
Ref. |l^ for a detailed discussion of its applicability to 
diamond surfaces. The use of the fully relaxed slab coor- 
dinates guarantees that the optical properties we calcu- 
late include all the subtle effects of surface- induced strain 
and appropriate geometry. 

To calculate the optical properties of the system, we 
generate the electronic level structure of the slab using a 
well known parameterized tight-binding approach with a 
sp^s* orbital basis.Ea This basis provides a good descrip- 
tion of the conduction band of cubic materials. This 
approximation has been applied to calculate, tb^ optical 
properties of a variety of silicon surfa(ie£,E3ll3 and re- 
cently to the (001) surface of diamond,E3 yielding good 
results. The parameters for CC interactions are taken 
to be the same as those of Ref. for the bulk, except 
for the on-site energy of the Pz orbitals of the surface 
atoms, Ep. This parameter is set to be 2.3 eV smaller 
than the corresponding bulk parameter. This change 
is assumed to be the likely result of additional orbital 
confinement at the surface, and as we will see below, 
it yields a level structure more attuned to experiments 
and other theoretical calculations. Moreover, the scaling 
factor of all tight-binding parameters for this particular 
surface was taken as (r/ro)^, where r is the bondlenght 
of any two first-neighbors atoms and To = 1.56 A, is 
the bondlenght in bulk diamond. These changes to the 
original bulk parameters provide an excellent description 
of the electronift-stcucture, as compared to experimental 
measurementspOflQ as we will show on Section III. 
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The optical properties of the surface region are deter- 
mined by its dielectric function. The imaginary part of 
the average slab polarizability is related to the transition 
probability between siab eigenstates induced by an ex- 
ternal radiation field.Ea Within a single-particle scheme, 
this relation is expressed by 



Im<iab(^) 



EElKc(k)|^ 



xSiEc(k)-E4-k)-hLu), (1) 



where p"^(k) is the matrix element of the a compo- 
nent of the momentum operator between valence (v) and 
conduction (c) states at the point k of the SBZ, 2c? is 
the slab thickness, m is the bare electronic mass, and 
A is the area of the surface unit cell. The real part 
of the average polarizability can be computed via the 
Kramers-Kronig relations. The surface dielectric tensor 
^surf(^) = 1 + 47ra"^f([jJ^^is then calculated from the 
average slab polarizability 
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Here abuik(w) = [ebuik('^) — l]/47r is the bulk polarizabil- 
ity, and dsurf is the depth of the surface region. Note that 
for cubic materials (C, Si, and Ge, for example) the bulk 
dielectrics-function is isotropic. The "tliree-layer model" 
of Drudctil and Mclntyre and Aspncf£3 adopted here is 
widely used in the analysis of optical data, and assumes 
that the system consists of three homogeneous regions: 
bulk, surface and vacuum, and the dielectric response is 
treated accordingly. 

The matrix elements of the momentum operator p"^(k) 
of Eq. (^, were obtained in terms of the atomic-like or- 
bital basis using the commutation relation between the 
Hamiltonian and position operator, p = i{m/h) [H,r]. 
Taking advantage of the orthogonality and localization 
of the orbitals, only the intra-atomic dipole matrix el- 
ements are retained. Then, only two additional pa- 
rameters to those of the tight-binding Hamiltonian were 
needed in order to reproduce the bulk dielectric func- 
tion. These parameters are the so-called intra-atomic sp 
and s*p dipoles, with best fitted values of 0.18 A and 
0.7 A, respectively. Notice, that these calculations-, ne- 
glect in principle excitonicEj and local field effects,t3 al- 
though the fitting parameter procedure compensates to 
some extent and yields very good agreement with bulk 
optical measurements. For a detailed description of the 
method the |Ceader is referred to the pioneering work of 
Selloni et alS^ and the review by Del SoIc.eZI 

In the above discussion, we have seen that the atomic 
structure of the surface region is intimately related to 
the dielectric response through its electronic structure, 
as given by Eq. (|^) . Experimentally, it is known that the 
surface dielectric function can be extracted by means of 
electronic and optical spectroscopies. Measurements of 
the reflectance anisotropy (RA) is one of these optical 
techniques which consists of measuring the relative re- 
flectance difference of two orthogonal light polarizations 



on the surface plane, x and y for example. Although 
the sample penetration of light is in general a few hun- 
dred times larger than the depth of the surface layer, 
the contribution from the bulk region to the RA spec- 
tra is canceled since the bulk optical properties of cubic 
materials are isotropic. Correspondingly, this technique 
is extremely sensitive to surface features and electronic 
properties due to reconstructions or adsorption events. 

Theoretically, the reflectivity is relatecLtp the dielec- 
tric function through the Fresnel formula]^ which must 
however, be modified due to the presence of the recon- 
structed surface regionETa This correction yields the 
following expression for the differential reflectance spec- 
trum isijien the light incidence is normal to the surface 
plane,! 



ARy 

Ro J 



Im 



ebuik(i:^) - 1 



(3) 



Here, a is one of the orthogonal directions on the surface 
plane, AR — R — Ro is the difference between the actual 
reflection coefficient R and the reflectivity Ro given by 
the Fresnel formula. 

The second experimental technique in which we arc in- 
terested is the electron energy loss spectroscopy (EELS). 
Here, an electron beam of a given low-energy and momen- 
tum is scattered by the sample. The electron beam in- 
duces polarizations on the surface region so that the elec- 
trons lose some of this energy before being scattered into 
the detector. The process can -be described well in terms 
of a dipolar scattering theoryj^j and provides a suitable 
description of vibrational modes of surface atoms and 
molecules, as well as electron transitions on the surface 
region. In the present work, all of the electronic transi- 
tions in the surface region are due to the reconstruction 
of the surface and not to adsorbates, although the work 
could be generalized to include various adsorbate species 
as well. 

The electron scattering probability P(q||,a;) for an 
electron that loses a quantum of energy hw, and transfers 
a momeniiim ?iq|| in the direction of the surface plane is 
given hyvB 



1 k' 



111 



(eao7r)2 cos (^i k \qj + q'l^l'^ 



Im5(q|h'^), (4) 



where k and k' are the wave- vectors of the incident and 
scattered electrons, ipi is the angle of incidence and hq± = 
h{kz — k'^) is the momentum transfer in the direction 
perpendicular to the surface plane. The above relation 
holds when the energy loss and momentum transfer to the 
medium are small. Assuming that the scattering occurs 
in the yz plane, the loss function is deflned by 



Im g{qy, uj) = Im 



1 + eoff(gy,t^) 



(5) 



where ece{qy,uj) is the nonlocal effective dielectric func- 
tion of the system. In the limit qzdgmf ^ 1, when the 
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momentum transfer to the medium in the perpendicu- 
lar direction to the surface plane is small, this effective 
dielectric function becomes 

ees{qy,U!) « ebuik(w) 

+ lydsurt [eZfi^) - ebulk('^)/<urf (^)] • (6) 

This theory has been applied successfully to explain the 
experimental EELS spectra of tlifi_2xl and 7x7 recon- 
structions of the Si(lll) surface.EjO 

In the following section we use our calculated surface 
dielectric function to explain the main features of RA 
and EELS spectra of the C(lll)-2xl surface. 



III. RESULTS AND DISCUSSION 

A. Surface Band Structure 

The electronic band structure of the C(lll)-2xl sur- 
face is presented in Fig. ^. The electronic structure is 
shown along the main symmetry directions of the irre- 
ducible SBZ, from F to J (x direction), from J to K 
{y direction) and from K to F (diagonal). The states 
associated to the surface reconstruction are represented 
by stars, while dots correspond to the projected bulk 
states. The top of the bulk valence band is set at eV, 
and the calculated Fermi level {Ef) is at about 1.5 cV 
(not indicated in Fig. |^), and coincident with the nearly- 
degenerate and flat dispersion states along JK. This 
result is in excellent agreement pJijzitlL the reported ex- 
perimental value of 1.5±0.2 eV.&ElQl3 In fact, the cal- 
culated results presented here are in excelleiit.|^neral 
agreement with experimental measurementsDLiij'lj and 
also compare welkYith those calculated previously us-, 
ing first-principlesErt2l and parameterized tight-bindingtll 
approaches. 

The calculated surface band structure of Fig. ^ shows a 
large gap of about 5.5 eV between the occupied and unoc- 
cupied surface states at the F point. The occupied surface 
states liying within the bulk valence band show a disper- 
sion of ~ 2.4 eV along the FJ direction (the Pandey chain 
axis direction), and have mainly a pz character. The be- 
havior of these states is similar to op.e observed experi- 
mentallw by Himpsel and co-workersH and Pate and co- 
workers, where a nearly flat filled surface band is found 
from F to about 0.5FJ where there is a minimum, and 
then rapidly disperses upward while approaching the J 
point. A similar behavior is found for these surface states 
along the FK direction. Above the Fermi level there are 
two bands of unoccupied states near 4.5 and 5.5 eV at the 
F point. These states have a strong s and p^ component 
corresponding to the dangling bonds of the surface chain 
atoms and the backbonds with the second and third layer 
atoms. The states at 4.5 eV show a nearly flat band in 
the first half of the FJ and FK directions, at about the 
halfway point in both directions the band has a maximum 
and then rapidly disperses downward approaching the J 



and K points. In the direction perpendicular to the chain, 
JK, these empty surface states and the occupied surface 
band become nearly degenerate and show little dispersion 
of less than 0.1 eV. Notice that these two states cross the 
Fermi level halfway through the JK direction. The strik- 
ing difference in dispersion of the surface bands along 
the two main directions is a reflection of the Pandey-like 
chains formed on the surface. The chains along the FJ 
(or x) direction allow for nearly free-like electronic mo- 
tion (although the dispersion is not parabolic), while the 
nearly vanishing overlap between the chains along JK (or 
y) direction yields nearly flat surface bands and reduces 
electronic hopping across zigzag chains. 

Near ~5.5 eV at the F point begins a band of unoc- 
cupied surface states mainly due to the surface dangling- 
bonds, and part to second layer backbonds, correspond- 
ing approximately to those calculateAJpy Vanderbilt and 
LouieO and Alfonso and co- workers, t2l and likely to be 
those observed by Kubiak and Kolasinski.D These states 
show a dispersion of about 1 eV with a minimum at about 
halfway the FJ direction, where they anticross the sur- 
face band associated to the backbond states described 
above, and produce a hardly noticeable splitting at the 
crossing. On the other hand, these states show less dis- 
persion ('^0.5 eV) along the FK direction and never cross 
the empty surface dangling bonds band. These states 
have also been observed experimentally by Kubiak and 
KolasinskiQ with a weak intensity at an energy of about 
5.8 eV from the top of the valence band, for both the FJ 
and FK directions. Finally, some localized resonance- like 
occupied states are also found at —2.5 eV and at —4.8 eV 
near the F point (shown as stars within the valence band 
in Fig. 1^). The former states are mainly due to the sub- 
surface chains with a strong p^, component, and the latter 
one have first- and second-layer backbond characteristics. 
The states at —2.5 eV, are similar to those reported by 
Vanderbilt and Louie E |-|r-, 

Experimentaljjjihotoemission results show occupiecH'B 
and unoccupiedQu surface states with a gap of nearly 
5.1 eV at the F point. The Fermi level is reported at 
about 1.5 eV above the top of the valence band, simi- 
lar to our findings, while the dispersion of the observed 
surface states is in very good agreement with those calcr^ 
lated here, and those reported by Vanderbilt and Louie.El 
Moreover, since these surface states are detectable only 
for p-polarization, they show a strong s and p^ charac- 
ter, in agreement with our results. A resonance unoccu- 
pied state has also been observed at about 6 eV from the 
top of the valence band at the F point .0 This resonance 
state is much weaker than the surface states lying in the 
fundamental gap, and it has not been possible to fully 
investigate its orbital character, although it could be the 
higher-energy surface state we find. 

While the calculated surface bands-ijt.Fig. ^ compare 
well |J4pt|h the experimentaL results, &BErlII the calculated 
LDAOIIj and tight-bindingliiJ results are rigidly shifted by 
-f 1 and +2 eV, respectively. The discrepancies among the 
different approaches could perhaps be partly attributed 
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to many-body effects. For example, when the exchange 
correlation effects are considered, the surface band is 
shifted Jtowards the top of the projected bulk valence 
bandE'El Moreover, when dynamical effects are taken into 
account within the GW approximation, the surface band 
moves into the projected bulk valence band in the vicin- 
ity of thje,r point also in agreement with experimental 
resultsBB However, there could be other sources of er- 
ror when we compare directly with experimental results, 
including the precise experipnental location of the Fermi 
level, as pointed out before.El 

The total electronic density of states (DOS) of the slab, 
and the projected density of states of the first two layers 
are shown in Fig. |. The DOS was calculated taking an 
average over 4900 points distributed homogeneously in 
the irreducible SBZ. We observe within the fundamental 
bulk gap, between and 5.5 eV, a non-zero density of 
states coming mainly from the dangling bonds associated 
with first layer atoms. This continuum of surface states 
is responsible for the metallic-like behavior of the surface 
around the Fermi level, as we will discuss in detail in 
the next section. Notice that the DOS from to 4 eV 
is nearly a constant, as expected for a 2D free-electron 
system. From Fig. ^ is clear that the peak at about 4.5 eV 
with a strong pz component has its origin in the dangling 
bond and the first- and second-layer backbonds. The 
resonance states at about —2.5 and 5.5 eV are associated 
with the second layer chains and the backbonds between 
the first layer and second layer atoms, respectively. The 
pronounced peaks of the projected DOS in the first layer, 
at about —1.5 and 4.5 eV are due to the lack of dispersion 
of the surface bands on the first half of the FJ and FK 
directions of the SBZ (see Fig. ^). Here, we observe that 
the states in the bulk gap are mainly localized in the first 
two layers, as one would expect, with decreasing intensity 
into the slab. 

Finally, before addressing the optical consequences of 
this level structure, we should comment on our choice of 
parameters. The excellent agreement with experiments 
and ab initio electronic calculations has been greatly en- 
hanced by our use of the different Ep^ parameter at the 
surface atoms, as mentioned above, as well as to the fully- 
relaxed atomic positions for the reconstructed surface of 
Ref. ^. Indeed, use of the bulk Ep^ parameters for all 
surface atoms yields a level structure (riot shown) very 
similar to that of Davidson and Pickett liil In that case, 
we obtain Ej « 3.5 eV above the valence-band top, while 
the filled surface dangling bond state remains w 2.7 eV 
below Ef (but now above the valence band). Similarly, 
the gap between surface state and conduction band bot- 
tom along the JK direction is only « 2.5 eV, rather than 
the 4 eV gap shown in Fig. 2. This full set of results vali- 
dates the choice of the physical parameter Ep^ at the sur- 
face. Although a detailed fit to the experimental results 
was not performed, it is clear (as one would expect on 
general physical grounds) that the orbital localization at 
the surface affects the diagonal tight-binding parameters. 
A full ab initio determination of the various optical pa- 



rameters, both in the bulk and near the surface, together 
with the fully relaxed level structure will be obviously de- 
sirable. We are currently carrying out such project and 
our results will be presented elsewhere. 



B. Surface Dielectric Properties 

The imaginary part of the average polarizability of the 
slab, Eq. (|^), was calculated using 4900 points distributed 
homogeneously on the irreducible SBZ. The large number 
of points needed is due to the small (large) size of the 
surface unit cell in real (reciprocal) space and to the large 
sections of the SBZ with flat joint density of states. The 
average over this large number of points is necessary to 
give full and reliable convergence of the optical properties 
for this particular surface. Electron transitions up to 
20 eV were taken into account, so that after the Kramers- 
Kronig transform the calculated real part is accurate up 
to about 10 eV. 

In Fig. m we present the real and imaginary parts 
of the surface dielectric tensor esurf(w) calculated from 
Eq. (H). The thickness of the surface region used was 
dsurf = 2.5 A, which approximately corresponds to two 
monolayers (other choices of dsurf do not change qual- 
itatively our results for energies in the bulk gap). The 
response for light polarized along the chains (x axis, Cgui-f ) 
corresponds to the solid lines, while for light polarized in 
the y direction (egu^f) shown by dotted lines, and the 
dashed lines correspond to the direction perpendicular to 
the surface plane (Cguif). The imaginary part of esuif(<-^) 
along X shows a strong peak at about 0.1 eV that is 100 
times more intense than the rest of the structure shown 
in Fig. ^. This peak at low-energy is a reflection of the 
metallic-like character of the surface along the chains. 
Then, from 2 to 5.5 eV the dielectric function is nearly 
constant up to the point when electron transitions be- 
tween bulk states become important. 

The following discussion about the origin of the main 
electron transitions of the surface dielectric function can 
be seen clearly in the lower four panels of Fig. where 
the reflectance anisotropy spectrum has been decom- 
posed into the different contributions, S-S, S-B, B-B, 
and B-S. For light polarized in both directions, x and y, 
the dielectric response is dominated by transitions among 
surface states (S-S) up to ~ 4 eV. From about 4 eV the 
contribution of the transitions from surface to bulk states 
{S-B) and from bulk to surface states (B-S) becomes im- 
portant. Note that the gap between the occupied flat 
band of surface states along the JK direction on the SBZ 
(see Fig. |^), and the bottom of the bulk conduction band 
is about 4 eV. Likewise for the gap between valence bulk 
states and the unoccupied surface band at F point. The 
high density of surface states above and below the Fermi 
level along FJ results in large S-B and B-S contributions 
to the Imesui.f(w) along the y direction (perpendicular to 
the chains; shown dotted in Fig. 2). Only for this per- 
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pendicular direction to the chains, the Imesurf(w) shows 
a intense peak centered at about 6 eV due to S-S tran- 
sitions. The transitions between bulk states (B-B) be- 
come important from about 5.5 eV onwards, where the 
response to x and y polarizations is very similar, as one 
expects for cubic semiconductors (notice also the scale 
change as the traces are much weaker). For light po- 
larized perpendicular to the surface plane {z direction), 
the eg^j.£ shows also a peak around 6 eV mainly due to 
the first- and second-layer backbond states (figure not 
shown) . 

In the rest of this section we will discuss the reflectance 
anisotropy and electron energy loss spectra obtained us- 
ing the calculated surface dielectric function. 

1. Reflectance Anisotropy 

In Fig. H, the top panel shows the differential re- 
flectance anisotropy spectra for light at normal incidence, 

(^) ^ (^) ' ^^^^^^^^^'^ according Eq. and la- 
belled TOTAL. This has been decomposed in its differ- 
ent contributions, where the response to light polarized 
along X (chain axis) corresponds to the solid line, while 
the dotted line corresponds to light polarized along y 
(perpendicular to the chain). From the figure, it is clear 
that the spectrum shows a large surface anisotropic op- 
tical response in a large range of photon energies. While 
for x-polarized light the spectrum shows mainly one peak 
at low energies, the y-polarization spectrum shows a rich 
structure for all energies inside the bulk optical gap. 

The intense peak at ~0.1 eV corresponding to x- 
polarized light is totally determined by transitions be- 
tween surface states. As we have explained above, this 
peak is related with the metallic-like behavior of the 
surface along the chain axis. At about 6 eV there are 
also some S-S transitions of weaker intensity for x polar- 
ization, and associated with the resonance states in the 
conduction band. The rest of the a;-polarized spectrum 
shows a very small contribution from S-B and B-S tran- 
sitions, compared to the response for y-polarized light. 
In fact, the response to light polarized perpendicular to 
the chain axis (y direction) shows much more structure 
in a larger energy region within the bulk optical gap. Up 
to ~4 eV the spectrum is only dominated by S-S transi- 
tions. At 4 eV the contribution from S-B and B-S tran- 
sitions starts and is reflected in the TOTAL differential 
spectrum by a shoulder. As mentioned above, 4 eV cor- 
responds to the gap between the flat band surface states 
around Ef along JK and the bottom of the conduction 
band, as well as to the energy difference between valence 
bulk states and the unoccupied surface band beginning at 
4.5 eV at the F point. Then, the intensity enhancement 
of the S-B contribution starting from ^ 6.5 eV corre- 
sponds to an increase of the density of the conduction- 
band states. In all cases, the B-B contributions to the 
TOTAL differential reffectance spectra in this range are 



insignificant, since both polarizations yield nearly iden- 
tical contributions. 

It is important to notice that this kind of deconvolu- 
tion of the spectrum helps one gain useful insights into 
the nature of the various transitions. As we have pointed 
out, the S-B contribution starts at some determined en- 
ergy (« 4 eV), as this gap is related with the conduction 
band and surface states located around Ef. This part of 
the spectrum gives then unambiguous information on the 
position of the Fermi level with respect to the bulk band 
structure, and therefore the energy at which the filled 
surface states are. Notice that one important advantage 
of this optical spectroscopy is the high precision in mea- 
suring the energy at which the electronic transitions oc- 
cur. The present results could be important for a future 
comparison with reflectance anisotropy measurements in 
order to better determine the electronic structure associ- 
ated to this particular reconstruction of the surface. We 
hope this motivates additional experiments. 

2. EELS 

The calculated scattering probability of an EELS ex- 
periment, using Eq. (^), is shown in Fig. ^. The primary 
energy of the electron beam was taken equal to 80 eV, 
with a normal incidence geometry. The spectrum corre- 
sponding to an electron beam polarized along the chain 
(solid line) is very different at low energies (less than 
1 eV) than the results for a polarization perpendicular to 
the chains (dotted line) . The intensity of the x-polarized 
reflected beam is a few thousand times larger that the 
beam for y direction polarization. At higher energies, 
from about 4 eV onwards, the two spectra show similar 
amplitude and behavior. The inset shows the scattering 
probability for energy loss from 2 to 8 eV, where the in- 
tensity has been augmented 500 times. Here, a feature 
starts at about 4 eV. As we discussed above, these struc- 
tures are related to the contribution to the surface di- 
electric response from S-B and B-S transitions, while the 
broad peak at about 6 eV is produced by S-S transitions. 
Notice that the high energy-loss features are strongly re- 
duced by the decaying prcfactor ^/q\\ appearing in the 
definition of P, Eq. (4). r. 

The EELS experiments reported by PepperQ showed a 
broad structure centered near 2.1 eV (and width of about 
1.7 eV). The primary energy of the normal incident elec- 
tron beam was Eo = 80 eV. The main spectrum reported 
by Pepper was obtained by subtracting the spectra mea- 
sured for the clean and hydrogenated surfaces, in order 
to reduce the effects of a strong elastic peak and to en- 
hance the signal due to the reconstruction. The spectra 
of the clean and hydrogenated surfaces were obtained 
from an average over the SBZ. In the difference spec- 
trum a minimum gap of about 1 eV was observed and 
identified with the effective gap between surface states 
at the point J. The energy resolution of the system is 
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estimated at 0.63 cV, from the width of the elastic peak 
remnant. Since the energy resolution is not optimal in 
this experiment, it is difficult to make a direct compari- 
son with theoretical calculations and other experiments. 
Moreover, the energy loss measured by this kind of spec- 
troscopy is generally smaller than the energy difference 
between occupied and unoccupied states of the system 
on its ground state. Therefore, it is possible that the ob- 
served broad feature at 2.1 cV is related to the S'-S' tran- 
sitions integrated over the SBZ, and expected to have 
an enhanced joint DOS at « 6 eV. The overall result- 
ing feature would perhaps be a combination of excitonic 
downshift and the high-energy {l/q\\) suppression factor. 
It is clear, however, that a better-resolution and more de- 
tailed EELS study on this system will be highly desirable. 
We will be glad to provide details of our electronic struc- 
ture and surface dielectric function results to interested 
experimental groups. 

IV. CONCLUSIONS 

We have investigated the optical response of the 
C(lll)-2x 1 surface based on a sp^s* parameterized tight- 
binding approach. The dielectric function of the surface 
region was calculated and a large anisotropy was found. 
This anisotropy of the optical response is a direct con- 
sequence of the surface reconstruction. The dielectric 
response of the surface was analyzed in terms of the S-S, 
S-B, B-B and B-S transitions, and important features 
corresponding to each type of transitions were found. 
The reflectance anisotropy and electron energy loss spec- 
tra were calculated in order to provide direct compari- 
son with experiments. We can conclude that these op- 
tical spectroscopies combined with theoretical studies, 
can help one elucidate the controversial surface electronic 
structure, and therefore, the structural and electronic 
level reconstruction of this important surface. 
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FIG. 4. The (a) imaginary part and (b) real part of the 
surface dielectric response. Solid lines correspond to light 

polarized along the chains (x axis) tg^j.^ , while dotted lines 
correspond to light polarized perpendicular to the chain {y 
axis), Cgurf The dashed line corresponds to light polarized 
perpendicular to the surface, e^^^f . 

FIG. 5. TOTAL Differential Reflectance spectrum divided 
into its S-S, S-B, B-B, and B-S components. Solid lines in 
bottom four panels correspond to light polarized along the 
chain {x axis), while the dotted lines correspond to light po- 
larized perpendicular to the chains {y axis). 

FIG. 6. Scattering probability function for electron energy 
loss, Eq. (4). Solid lines correspond to light polarized along 
the chain (x axis), while dotted lines correspond to light po- 
larized perpendicular to the chains {y axis). The inset has 
been amplified 500 times. 



FIG. 1. Atomic model of the C(lll)-2xl surface, (a) Top 
view with the three uppermost layers; dashed line corresponds 
to the surface unit cell, (b) Side view with the six upper- 
most layers. The first-layer atoms forming Pandoy-liko chains 
shown in black, (c) Surface Brillouin Zone is shown; shadowed 
area corresponds to its irreducible part. The main symmetry 
points are indicated. 



FIG. 2. Surface electronic structure along the main sym- 
metry directions of the surface unit cell. Dots correspond to 

the projected bulk states, while stars represent surface states. 
Resonance states embedded in projected bulk bands are also 
represented by stars. 



FIG. 3. Calculated density of states for (a) TOTAL, (b) 
projected on first layer 1st LAYER, and (c) projected on 2nd 
LAYER (solid) and 3rd LAYER (dotted fine). Different ver- 
tical scales used in each panel. 
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Fig. 1 "Anisotric optical response of diamond (111)-2x1 surface" 
C. Noguez and S.E. Ulloa 
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Fig. 2 "Anisotropic optical response of the diamond (111)-2x1 surface" 
C. Noguez and S.E. Ulloa 
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Fig. 3 "Anisotropic optical response of diamond (111)-2x1 surface" 
C. Noguez and S.E. Ulloa 
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Fig. 5 "Anisotropic optical response of diamond (111)-2x1 surface" 
C. Noguezand S.E. Ulloa 




Fig. 6 "Anisotropic optical response of diamond (111)-2x1 surface" 
C. Noguezand S.E. Ulloa 



